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The  dislocation  distribution  as  a function  of  depth  was  investi- 
gated on  fatigued  and  unidirectionally  strained  specimens  of  4l"50 
steel,  2024  aluminum  as  well  as  on  single  crystals  of  aluminum 
and  silicon.  Using  X-ray  line  broadening  together  with  the  Berg- 
Barrett  technique  it  was  shown  that  the  excess  dislocation  den- 
sity was  much  larger  in  the  surface  layer  to  a depth  of  about 
100  micrometres  than  in  the  bulk.  From  the  measurements  of  the 
initial  permeability  at  high  frequencies  (10  kilohertz)  it  would 
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also  be  inferred  that  the  dislocation  density  was  larger  in  the 
surface  layer.  Because  acoustic  emissions  were  found  at  very  low 
stresses  relative  to  the  maximum  stress  during  fatigue,  it 
appears  that  this  method  cannot  be  used  on  ferromagnetic  mate- 
rials to  detect  fatigue  damage 
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ABSTRACT 


The  dislocation  distribution  as  a function  of 
depth  was  investigated  on  fatigued  and  unidi- 
rectionally  strained  specimens  of  4130  steel, 

2024  aluminum  as  well  as  on  single  crystals 
of  aluminum  and  silicon.  Using  X-ray  line 
broadening  together  with  the  Berg-Barrett 
technique  it  was  shown  that  the  excess  dislo- 
cation density  was  much  larger  in  the  surface 
layer  to  a depth  of  about  100  micrometres 
than  in  the  bulk.  From  the  measurements  of 
the  initial  permeability  at  high  frequencies 
(10  kilohertz)  it  would  also  be  inferred  that 
the  dislocation  density  was  larger  in  the 
surface  layer.  Because  acoustic  emissions 
were  found  at  verv  low  stresses  relative  to 
the  maximum  stress  during  fatigue,  it  appears 
that  this  method  cannot  be  used  on  ferromag- 
netic materials  to  detect  fatigue  damage. 

INTRODUCTION 

It  has  been  shown  previously1  that  during  uniaxial  fatiguing 
of  metals  the  work  hardening  of  the  surface  layer  increases  more 
rapidly  than  that  in  the  bulk.  For  high-cycle  fatigue  after  the 
first  1/2  cycle  (R  = 1)*  the  work  hardening  was  confined  entirely 
to  the  surface  layer.  From  a systematic  investigation6  on  alum- 
inum 2014-T6,  titanium  (6  Al/4v)  and  a 4130  steel  it  was  observed 
that  a propagating  fatigue  crack  was  formed  whenever  the  surface 
work  hardening,  as  characterized  by  the  surface  layer  stress,  Ts* 
reached  a critical  value.  This  critical  value  was  independent  of 
the  stress  amplitude,  prior  fatigue  stress  history,  and  environ- 
ment. These  parameters  affected  the  rate  of  surface  work  harden- 
ing but  not  the  critical  value  for  the  fatigue  crack  formation. 
The  observation  on  crack  formation  led  directly  to  a mechanism 
for  fatigue  damage.  It  was  proposed6  that  since  the  surface 
layer  acted  to  oppose  the  motion  of  dislocations  it  could  provide 
a barrier  to  support  a piled  up  array  or  accumulation  of  disloca- 
tions of  like  sign.  When  the  barrier  is  of  sufficient  strength 
to  support  a critical  number  of  these  dislocations,  fracture  will 
occur  when,  locally,  the  stress  fields  associated  with  the  accum- 
ulation of  the  dislocations  exceed  the  fracture  strength  of  the 
material.  Accordingly,  fatigue  damage  may  be  described  in  terms 
of  the  surface  layer 


1 Superscripts  refer  to  similarly  numbered  entries  in  the  Techni- 
cal References  at  the  end  of  the  text. 

♦Definitions  of  abbreviations  used  are  on  page  i. 
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where  N0  is  the  number  of  cycles  to  form  a crack  and  ts  is  the 
work  hardening  of  the  surface  layer. 

However,  since3 


where  a is  a constant  of  the  order  of  1/5#  G is  the  shear 
modulus,  b is  the  burgers  vector  and  p is  the  dislocation  den- 
sity. Since  Taira,  et  al,  have  shown  that  the  ratio  of  the 
excess  dislocation  of  one  sign  , D,  to  the  total  dislocation 
density,  p,  is  a constant,  then 


From  equations  (3)  and  (4)  it  appears  that  the  extent  of  fatigue 
damage  may  be  determined  by  measuring  the  dislocation  density  in 
the  surface  layer.  Therefore,  to  this  end,  three  methods  for  the 
nondestructive  evaluation  of  the  dislocation  density  in  the  sur- 
face layer  were  investigated: 

• X-ray  line  broadening. 

• Magnetic  characteristics. 

• Acoustic  emission. 

For  the  X-ray  portion  of  this  investigation,  silicon  and 
aluminum  were  selected  as  model  materials  to  study  the  effect  cf 
plastic  deformation  and  fatigue  cycling  on  surface  layers  and 
bulk  material.  Silicon  may  be  regarded  as  a representative  of  a 
class  of  materials  having  low  stacking  fault  energy, 6 9 whereas 
aluminum  is  a typical  representative  of  materials  with  high 
stacking  fault  energy.  The  stress-strain  curves  of  silicon 
and  other  diamond  structure  materials,  when  taken  at  elevated 
temperatures,  show  a pronounced  analogy  to  fee  metals.9  However, 
at  low  deformation  temperatures,  silicon  is  extremely  brittle  and 
exhibits  the  deformation  characteristics  of  a typical  ceramic 
material.  Thus,  if  dislocations  are  introduced  into  crystals  at 
elevated  temperatures,  the  dislocation  configuration  remains 
"frozen- in"  when  the  temperature  is  lowered  to  room  temperature 
since,  at  ambient  temperatures,  the  dislocations  are  virtually 
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immobile.  Consequently,  in  studying  the  deformation  response  as 
a function  of  depth  from  the  surface,  removal  of  surface  layers 
by  chemical  polishing  is  not  likely  to  disturb  the  induced  dis- 
location configuration.  While  silicon  and  aluminum  single  crys- 
tals were  subjected  only  to  tensile  deformation,  the  aluminum 
2024  alloy  was  subjected  to  both  tensile  deformation  and  fatigue 
cycling. 

It  has  been  shown  in  a detailed  theoretical  study1 0 that 
dislocations  affect  the  magnetic  properties  via  the  interactions 
of  the  stress  field  surrounding  the  dislocation  with  the  magnet- 
oelastic forces  within  each  crystalline  grain.  Although  the 
derivation  of  the  various  equations1 0 is  strictly  applicable  only 
at  or  near  magnetic  saturation,  it  is  indicated  that  dislocation 
strains  affected  other  portions  of  the  magnetization  curve  as 
well.  A more  simplified  derivation11  reproduced  in  appendix  A 
shows  that  whenever  magnetization  is  due  to  independent  domain 
rotation  then 


s Ti 


(5) 


where  is  the  initial  permeability.  Is  is  the  saturation 
mangetization,  Xs  is  the  magnetostrictive  coefficient,  and  Ti 
is  the  average  value  of  the  internal  stress.  From  equation  (2) 


Po  - 1 = 


Is 

XsabCps 


i/2 


(6a) 


and  since 


p/D  = C 


(6b) 


then 


u - 1 = s 

° XgabCGDs1/2 


(6c) 


Over  and  above  the  parameters  usually  associated  with  the 
magnetization  of  a ferromagnetic  metal,  the  phenomenon  of 
Barkhausen  noise1  has  recently  been  fpund  to  be  associated  with 
both  residual  stress1,3  and  grain  size1  in  steel.  The  Barkhausen 
effect  has  been  established  as  being  due  to  the  relatively^ dis- 
continuous nature  of  the  mangetization  of  a ferrous  metal. 
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In  certain  regions  of  the  magnetization  curve,  favorably  oriented 
domains  (i.e.,  regions  of  uniform  magnetization)  grow  in  volume 
at  the  expense  of  those  least  favorably  oriented.  This  growth  in 
volume  does  not  take  place  continuously  (assuming  a continuously 
changing  magnetizing  field)  but  executes  many  small,  rapid  jumps. 
The  resultant  flux  changes  will  generate  many  transient,  voltage 
spikes  in  a suitably  disposed  pickup  coil.  Hence,  observation  of 
these  voltage  spikes  should  yield  some  indication  of  grain  size 
and  dislocation  density.  In  summation,  there  is  plentiful  theo- 
retical and  experimental  evidence  that  suggests  that  measurement 
of  the  ferromagnetic  properties  of  steel  near  the  surface  should 
yield  a measure  of  the  surface  dislocation  density. 

Measurement  of  any  ferromagnetic  property  at  the  surface  of 
a steel  specimen  must  take  due  account  of  the  "skin  effect, " 
i.e.,  the  lack  of  penetration  of  a varying  magnetic  field  below 
the  surface  of  an  electrical  conductor.  This  propensity  for 
time-varying  fields  to  decrease  with  depth  in  conductors  is  due 
to  the  "eddy  currents"  created  by  the  electromotive  force  induced 
by  the  varying  magnetic  fields.  These  currents  create  a secon- 
dary magnetic  field  that  interacts  with  the  primary  field  and 
causes  the  net  field  to  decrease  with  depth.  The  treatments  of 
this  phenomenon  most  suited  to  the  needs  at  hand  appear  in 
several  treatises.16  16  17  As  shown  in  appendix  B,  starting  with 
the  basic  Maxwell  field  equations,  we  find  that  the  magnetic 
field  inside  a planar  and  a cylindrical  conductor  decreases  with 
depth  in  accordance  with  the  following  relationships 


H = HQ  e“x/s  cos  |2rr  ft  - gj  (7) 


where  x is  the  depth  below  the  surface  of  a semi-infinite  planar 
conductor  and. 


H _ h ber2f2Sr/d)  + bei2f2&r/d) 
ber2  9 + bei2  6 


(8) 


For  a cylindrical  conductor  of  diameter  d,  r is  the  distance  from 
the  center  of  the  cylinder  to  the  point  in  question,  and  ber  and 
bei  are  Bessel  functions.  in  either  case,  the  magnetic  field  at 
the  surface  is  parallel  to  the  surface  and  to  the  cylinder  axis 
and  is  assumed  to  vary  sinusoidally  at  a frequency  f.  Other  sym- 
bols are  as  defined  in  appendix  B.  The  planar  solution  yields 
some  important,  easily  interpretable  consequences  that  carry  over 
to  the  less  interpretable  cylindrical  solution.  it  can  be  shown 
that  the  magnitude  of  H decreases  to  1/e  of  its  value  at  the  sur- 
face when  the  depth  x is  equal  to  s which  is  defined  by 
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s = 


(9) 


(TT^Of  )-1/2 


The  phase  of  the  field  referred  to  the  surface  shi  ts  with  depth 
so  that  at  one  skin  depth  the  field  lags  the  surface  field  by  1 
radian.  The  implication  of  this  is  that  the  effective  depth  for 
determining  the  electromotive  force  produced  by  flux  changes  in 
a specimen  is  somewhat  less  than  the  skin  depth  determined  by 
equation  (9)-  Finally,  for  materials  of  given  magnetic  proper- 
ties, i.e.,  permeability  p and  conductivity  a,  the  skin  depth 
varies  inversely  as  J~E.  Hence,  the  higher  the  frequency,  the 
more  indicative  is  the  response  to  the  value  of  d and  of  o at  the 
surface,  and  measurements  at  several  frequencies  should  yield  an 
indication  of  differences  between  volume  and  surface  microstruc- 
ture and  residual  stress. 

A detailed  theory  of  the  effects  of  eddy  currents  has  been 
workedg out  largely  in  connection  with  their  use  as  an  inspection 
tool.1  '1?  Confining  our  attention  initially  to  the  geometry  of 
cylindrical  specimens,  we  find  that  the  effects  of  eddy  currents 
are  most  conveniently  expressed  by  introducing  the  concept  of 
effective  permeability.  It  enters  in  the  following  way:  assume 

a cylinder  of  diameter  d,  permeability  M,  and  conductivity  o. 

If  the  cylinder  is  surrounded  by  a time-periodic  but  spatially 
unvarying  magnetic  field  H,  then  the  electromotive  force  induced 
in  a coil  of  n turns  of  diameter  dc  will  depend  on  the  applied 
field  H,  the  permeability  of  the  cylinder,  and  the  magnitude  and 
phase  of  the  eddy  currents  in  the  cylinder.  The  resultant 
electromotive  force  will  be  out  of  phase  with  the  applied  field 
and  in  very  general  terms  can  be  shown  to  be  given  by 


^1-T!  + iWeff>r^ 


(10a) 


for  the  component  that  would  normally  be  present  in  the  absence 
of  eddy  currents  and  a quadrature  component  Ei,  where 


where 


(10b) 


(10c) 


^eff.r  is  the  real  component  of  the  effective  permeability, 
^eff.i  is  the  immaginary  component,  and  r)  = (d/do)2  is  the  so- 
called  fill  factor.  Note  that  the  magnitude  of  the  induced 
electromotive  force  E is  given  by 
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— 


E = ^Er2  + E i2J  /Ci 


(11) 


and  the  phase  angle  cp  between  applied  field  and  induced  emf  is 
given  by 


9 = tan“l 


(12) 


The  formalism  that  has  been  stated  by  equations  (10)  through 
(12)  takes  on  real  meaning  when  numerical  values  of  the  components 
of  ^leff  can  be  assigned.  These  have  been  calculated  from  basic 
considerations  arising  out  of  equation  (3)  as  outlined  in  appendix 
B.  The  resultant  numerical  values  are  embodied  in  figure  1.  7 


Figure  1 

Components  of  Effective  Permeability 
of  Various  Frequencies 
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The  key  parameters  for  entering  these  curves  are  the  fill  factor 
previously  defined  and  the  f/fg  ratio,  where  f is  the  frequency 
of  H and  fq  is  a characteristic  frequency  of  the  material  and 
geometry  given  by 


f 


9 


5066 

pod2 


Hz 


(13) 


where  the  units  of  o and  d are  metres/ohm-mm  and  cm,  respec- 
tively. The  values  of  ^eff  for  various  f/fg  have  also  been  tab- 
ulated,16 the  results  being  reproduced  herein  in  table  1. 


TABLE  1 

Effective  Permeability 
as  a Function  of  f/fg 


£/f<* 

^eff.r 

^ef f . i 

0 

1.000 

0.000 

1 

O.9798 

0.1216 

2 

0.  264 

0.2234 

5 

0.6992 

0.3689 

10 

0.4678 

0.34Q4 

20 

0.3180 

0.2657 

50 

0.2007 

0.1795 

100 

0.1416 

0.1313 

200 

0.1001 

0.0950 

400 

0.0707 

0.0682 

1000 

0.0447 

0.0437 

10,000 

0.0141 

0.0140 

For  considerations  closely  related  to  those  that  lead  to  equation 
(10),  the  inductance  and  resistance  of  a helical  coil,  concentric 
to  a long,  cylindrical  specimen,  are  given  by1 

L = L0  (l  - Tl  + Weff.r)  (l4a) 

for  the  inductance  L,  and 

R - Rq  = Lo  (^eff.lj  (l4b) 

for  the  resistance  R.  Note  that  L0  and  R0  are  the  coil  induc- 
tance and  resistance  when  devoid  of  a cylindrical  specimen,  and 
iu  =s  2TTf.  The  requires  values  of  t-^eff  maY  be  obtained  from 
figure  1. 
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From  the  foregoing  it  is  apparent  that  measurement  of  either 
the  amplitude  and  phase  of  the  induced  voltage  or  the  impedance 
of  a helical  coil  surrounding  a cylindrical  specimen  is  dependent 
on  the  electromagnetic  properties  of  the  specimen.  Furthermore, 
the  higher  the  frequency  of  measurement,  the  more  the  results 
should  depend  on  surface  properties,  since  both  magnetic  fields 
and  eddy  currents  fall  off  with  depth  in  a manner  approximately 
according  to  equation  (9).  Although  the  theoretical  equations 
for  skin  depth  are  simpler  for  planar  specimens,  the  relative 
ease  of  controlled  measurements  on  cylindrical  specimens  suggests 
their  use  initially.  As  for  the  other  parameters  of  interest, 
coercive  force,  Rayleigh  constant,  and  Barkausen  noise,  measure- 
ments were  not  accomplished  in  this  reporting  period,  so  that 
further  theoretical  discussion  will  be  deferred. 

The  possibility  of  detecting  fatigue  damage  by  acoustic 
emission  is  based  on  the  Kaiser  effect.  b Kaiser  reported  when 
specimens  of  steel,  aluminum,  copper,  pine  and  boxwood  were 
loaded  in  tension,  and  then  unloaded,  upon  reloading  no  acoustic 
emission  activity  was  detected  until  the  previous  maximum  stress 
was  exceeded.  Later,  Duke  and  Kline1  b reported  that  prestressed 
specimens  of  JO-30  brass  and  1100  aluminum  exhibited  acoustic 
emission  bursts  at  stresses  below  the  previous  maximum  stress 
whenever  the  surface  layer  was  removed.  This  observation  coupled 
with  the  observation  that  the  surface  layer  work  hardened  more 
than  the  interior  in  fatigue  as  well  as  in  uniaxial  deformation 
led  to  the  idea  that  a measure  of  fatigue  damage  could  be 
obtained  by  determining  the  increase  in  the  acoustic  emission 
threshold  stress  as  a function  of  the  number  of  fatigue  cycles. 
This  concept  appeared  feasible  provided  other  emission  sources, 
such  as  growth  of  magnetic  domains,  did  not  operate. 

EXPERIMENTAL  PROCEDURE 


X-RAY 


Single  crystals  of  silicon  and  aluminum  were  sliced  into 
thin  wafers.  Tensile  specimens  with  gage  dimensions  of  5 x 1.5 
mm  and  lengths  20  mm  (Si)  and  15  mm  (Al)  were  prepared  by  means 
of  a Servomet  spark-erosion  cutter.  The  silicon  specimens  had 
(112)  surface  orientation  and  [lio]  tensile  axis,  while  the 
aluminum  specimens  had  (100)  surface  orientation  and  [l00]  ten- 
sile axis.  Surface  damage  was  removed  by  lapping  with  600-grit 
silicon  carbide  and  by  chemical  polishing  (1  part  HF,  3 parts 
HNO3,  2 parts  CH3COOH)  for  the  silicon,  and  by  spark-erosion 
planing  followed  oy  electropolishing  (1  part  HNO3,  3 parts 
CH3COOH  for  aluminum).  At  least  70  Wn  was  removed  to  obtain 
damage-free  surfaces.  The  aluminum  specimens  were  annealed  at 
600°  C for  20  hours  in  a protective  argon  atmosphere.  All  spec- 
imens were  deformed  on  an  Instron  universal  tester.  Silicon 
specimens  were  deformed  at  a strain  rate  of  5*9  x 10”5  sec~l 
above  the  ductile-brittle  temperature  in  an  argon  atmosphere. 
Aluminum  specimens  were  tested  at  a strain  rate  of 
3.3  x 10-3  sec-1  below  0°  C and  maintained  at  0°  C during  subse- 
quent X-ray  analysis. 
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The  aluminum  2024  alloy  was  obtained  as  bar  stock  with  com- 
position (weight  percent)  3-99  Cu,  1.35  Mg,  0.50  Mn,  and  princi- 
pal impurity  content  of  0.29  Fe,  0.10  Si,  0.08  Zn,  0.02  Ti,  and 
0.02  Cr.  The  stock  was  fabricated  into  cylindrical  fatigue  spec- 
imens with  a gage  diameter  of  3 nun  and  a gage  length  of  6.35  mm. 
Specimens  were  heat  treated  in  a tilting  furnace  at  495°t5°  c for 
2 hours  in  an  argon  atmosphere,  quenched  in  ice  water  at  0°  C, 
and  aged  at  ambient  temperatures  for  48  hours.  After  heat  treat- 
ment the  specimens  were  polished  mechanically  with  polishing 
alumina,  using  0.05-dm  particles  for  the  final  step  in  the  pro- 
cess. Any  remaining  surface  damage  was  removed  by  electropolish- 
ing approximately  100  dm  from  the  gage  diameter  with  a 10$ 
perchloric,  90$  methanol  solution.  The  specimens  were  fatigued 
in  the  tension-compression  mode  on  a Tatnell-Krause  fatigue 
machine  with  constant  stress  amplitude  and  zero  mean  stress. 

The  applied  stress  amplitude  corresponded  to  one-half  of  the 
static  yield  stress,  and  the  number  of  cycles  (10^)  to  which  the 
alloy  was  subjected  corresponded  to  about  one-tenth  of  the 
expected  fatigue  life. 


The  X-ray  double  crystal  diffractometer  method  combined  with 
X-ray  reflection  topography  ( Berg-Barrett  Method,  hereafter 
referred  to  as  B-B  Method)  was  employed  as  the  principal  research 
tool  of  structural  analysis.  The  X-ray  rocking  curves  were  taken 
in  the  nx,  -ny)  arrangement  where  n refers  to  the  order  of  reflec- 
tion and  x ana  y to  the  first  crystal  and  test  crystal,  respec- 
tively. in  this  type  of  arrangement  the  normals  to  the  crystal 
surfaces  are  opposite  in  direction.  A dislocation-free  silicon 
crystal  in  the  (111)  orientation  was  used  as  the  first,  mono- 
chromating  crystal,  and  the  Cu  Kc*2  reflection  was  eliminated  with 
a knife  edge.  B-B  topographs  taken  at  angular  rotation  positions 
of  the  rocking  curves  allowed  for  correlation  between  the  config- 
uration of  the  defect  structure  and  detailed  features  of  the 
rocking  curves.60 “ The  lattice  defects  induced  by  fatigue 
cycling  of  the  A1  2024  alloy  were  studied  by  a double  crystal 
diffractometer  method  developed  for  application  to  polycrystal- 
line materials.  This  method  related  the  analysis  of  the  rocking 
curves  of  the  reflecting  grains  to  their  topography  and  mor- 
phology. 1 ' 3 - In  applying  this  method  each  reflecting  grain 
of  the  specimen  is  considered  to  function  independently  as  the 
second  crystal  of  a double  crystal  diffractometer,  witn  a 
reflected  beam  from  the  (111)  face  of  a silicon  crystal  as  the 
incident  beam  on  the  specimen,  and  carrying  out  angular  specimen 
settings  and  appropriate  film  shifts  between  the  settings  pro- 
duces pictorial  rocking  curves  for  each  reflecting  grain.  The 
rocking  curve  is  manifested  as  an  array  of  spots  of  rising  and 
falling  intensity  as  the  grain  is  rotated  through  its  angular 
range  of  reflection,  as  shown  by  figures  9»  10,  and  11  (see 
pages  15-17*  Berg-Barrett  topographs  taken  at  the  specimen 
surface  correlate  the  topography  and  morphology  of  the  grains 
to  the  corresponding  rocking  curves. 
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MAGNETIC  PERMEABILITY 

The  induced  voltage  measurements  were  performed  as  depicted 
in  figure  2. 


I,  4400  TURNS  OF  NO  30  WIRE  ON  9/16"  FORM,  2"  LENGTH,  15/16"  0 D 
12  - 100  TURNS  OF  NO  30  WIRE  ON  23/32"  FORM 


Figure  2 

Arrangement  for  Induced  Voltage  Measurement 
and  B/H  Curve  Trace 


The  alternating  magnetic  field  is  produced  by  driving  at  the 
desired  level  and  frequency.  Current  levels  in  are  estab 
lished  by  monitoring  the  potential  drop  across  the  current- 
sensing resistor  which  also  provides  the  horizontal  axis  of 
the  display.  The  voltage  induced  in  the  secondary  coil  Lp  is 
impressed  across  the  vertical  axis.  Vertical  axis  gain  is 
adjusted  so  that  both  horizontal  and  vertical  peak  deflections 
are  identical.  By  measuring  the  amplifier  gains  at  which  this 
occurs,  the  relative  amplitude  of  the  voltage  induced  in  L2  may 
be  determined.  The  resultant  elliptical  Lissajou  figure  seen  on 
the  oscilloscope  display  is  easily  converted  to  phase  angle  as 
shown  in  detail  in  appendix  C.  A Polaroid  photograph  of  the 
elliptical  trace  with  a superimposed  calibration  signal  preserves 
the  data  for  future  analysis.  A minor  modification  of  this 
arrangement  was  used  to  display  the  B/H  hysteresis  curve.  This 
was  accomplished  by  integrating  the  electromotive  force  induced 
on  the  secondary  L2  by  means  of  the  simple  RC  integrating  circuit 

MAT-76-74  10 


shown  in  figure  2.  The  mathematical  treatment  for  the  B/H  trace 
is  straightforward  and  is  given  in  appendix  D.  The  specimen  was 
demagnetized  prior  to  all  measurements  by  bringing  it  up  to 
saturation  at  low  frequency  (5  Hz)  and  then  gradually  decreasing 
the  amplitude  of  the  magnetic  field  until  it  vanished.  As  may 
be  seen  from  our  previous  discussion,  this  assured  measurement  of 
the  initial  permeability,  pQ.  Coil  impedance  measurements 
utilized  the  bridge  circuit  shown  in  figure  J>. 


SPECIMEN 


to  differential 
AMPLIFIER  TYPE  D 
plug-in  vertical 
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POWER 

AMPLIFIER 

TO  TEKTRONIX  536  OSC 
HORIZONAL  AMPLIFIER 


|.s  - DECADE  INDUCTORS,  SPRAGUE 
TYPES  850  WAI.2,3 

Rs-  DEKABOX,  MODEL  DB-265 

R - 5 fl,  PRECISION  NONINDUCTIVF 
WU  RESISTORS 


L,  SAME  AS  FIGURE  2 


Figure  3 

Low  Frequency  Impedance  Measurement  Arrangement 


The  coil  surrounding  the  specimen  was  the  same  as  Iq  in  figure 
2.  Current  levels  through  the  coil  were  maintained  in  the  linear 
region  of  the  B/H  characteristic.  The  standard  inductance  Ls 
and  resistor  Rg  are  adjusted  until  the  output  of  the  differential 
amplifier  as  displayed  on  the  oscilloscope  indicated  a null. 
Readings  of  the  Ls  and  Rg  settings  were  subsequently  corrected 
for  the  inevitable  small  resistive  component  in  Ls  and  inductive 
component  in  Rs. 
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ACOUSTIC  EMISSION 


The  acoustic  emission  measurements  were  made  at  the  Johns 
Hopkins  University  in  tensile  machinery  designed  for  quiet  oper- 
ation. The  load  was  applied  in  a continuous  manner  by  running 
water  from  a storage  tank  located  above  the  testing  machinery 
into  a container  at  the  end  of  the  lever  arm  (15:1  mechanical 
advantage).  A uniform  flow  rate  was  obtained  by  maintaining  a 
constant  pressure  head  in  the  storage  tank.  It  should  be  noted 
that  the  grip  configuration  allows  the  direct  placement  of  the 
transducers  on  the  specimen  without  contact  with  the  deforming 
region. 

The  monitoring  system  for  the  tests  consisted  of  a wide  band 
transducer,  preamplifier,  band-pass  filter,  and  counter.  A 
Parametrics  500  kHz  (critically  damped)  transducer  was  used.  The 
combined  amplifier  (Tektronix  1A7A),  preamp  and  counter  system 
used  has  the  capability  of  280-dB  gain  over  a frequency  range 
from  100  kHz  to  3 MHz,  ±3  db.  In  these  tests  the  signals  were 
limited  to  a frequency  range  from  75  to  200  KHz  by  a band-pass 
filter  (Kronhite  model  3202). 

SPECIMENS 

For  the  magnetic  portion  of  this  investigation  a 4130  steel 
specimen  having  a 12.5-mm  diameter  and  a 177.8-mm  gage  length 
was  used  for  the  tensile  and  fatigue  tests.  Bee  se  of  the 
adverse  ratio  of  length  to  diameter,  this  type  of  specimen  cannoL 
be  used  in  complete  reverse  axial  loading  without  buckling. 

These  specimens  were  normalized  at  815°  C in  1 hour.  The  speci- 
mens used  for  acoustic  emission  had  a 4-mm  diameter  and  a gage 
length  of  7 .6  nan.  They  were  austenized  at  87O0  C for  25  minutes, 
oil  quenched  and  tempered  at  5 60°  C for  30  minutes.  This  treat- 
ment gave  a 0.2#  yield  strength  of  156,000  psi  (IO76  MPa)  while 
the  normalized  specimens  had  a yield  strength  of  100,000  psi 
(689.5  MPa).  In  all  cases  after  heat  treatment  the  specimens 
were  mechanically  and  electrochemically  polished  to  remove  any 
decarburized  layer  that  may  have  been  formed. 

RESULTS 


X-RAY  ROCKING  CURVES 

A typical  halfwidth  value  of  the  rocking  curve  of  an  unde- 
formed silicon  crystal  was  about  15  seconds  of  arc.  After  10# 
tensile  deformation  at  65O0  C the  rocking-curve  profile  becomes 
multipeaked  as  shown  in  figure  4(a),  and  the  extent  of  the  half- 
width (3  pertaining  to  the  principal  peak  was  935  seconds  of  arc. 
Surface  layers  were  removed  by  chemical  polishing  in  steps  of  30, 
100,  175,  and  250  Pm.  Figure  4(b)  shows  the  rocking-curve  pro- 
file after  removal  of  surface  layers  totaling  250  Pm.  Both  the 
total  range  of  reflection  and  the  halfwidth  of  the  subpeaks 
decreased  drastically. 
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Figure  4 - X-ray  Rocking-Curve 
Profiles  of  Silicon  Single 
Crystal  Deformed  at  605°  C, 
e = 10$,  tensile  axis  [110]. 

(a)  Original  surface 
(b)  After  removal  of  25 0 dm 


The  dislocation  gradient  from  surface  to  bulk  was  calculated 
from  the  relation  given  by  Hirsch;  1 D = P^/9b  # where  D is  the 
number  of  excess  dislocations  of  one  sign,  0 is  the  measured 
halfwidth,  and  b is  the  magnitude  of  the  burger  vector.  The 
variation  of  dislocation  density  with  distance  from  the  surface 
is  shown  in  figure  5*  It  may  be  seen  that  after  removal  of  about 
100  dm  the  dislocation  density  was  virtually  that  of  the  bulk 
material . 


Figure  5 

Dislocation  Density  Gradient 
from  Surface  to  Bulk 
for  Silicon  Single  Crystal 


Tensile  deformation  of  aluminum  to  10$  strain  at  0°  C also 
resulted  in  a multipeaked  rocking  curve  as  shown  in  figure  6(a). 
The  lattice  misalignment  expressed  in  terms  of  the  0 value 
increased  fourfold  compared  to  the  annealed  crystal.  Figure 
6(b)  shows  the  rocking-curve  profile  after  surface  removal  of 
100  dm.  It  may  be  observed  that  the  total  lattice  misalignment 
as  well  as  the  halfwidth  values  decreased.  The  dependence  of 
the  ratio  of  dislocation  density  of  surface  layers  to  bulk  den- 
sity as  a function  of  depth  distance  is  shown  in  figure  J. 
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Figure  6 

X-ray  Rocking-Curve  Profiles 
of  Aluminum  Single  Crystal 
Deformed  at  0°  C,  e=  10$, 
tensile  axis  [100]. 

(a)  Original  surface 
(b)  After  removal  of  100  dm 


Figure  7 

Dislocation  Density  Gradient 
from  Surface  to  Bulk 
for  Aluminum  Single  Crystal 


Focusing  attention  on  the  B-B  topographs,  figure  8,  one  may 
note  that  the  average  size  of  the  substructural  reflecting 
domains  pertaining  to  the  surface  layer  was  smaller  than  that 
obtained  after  removal  of  a layer  of  100  dm. 


12  34 
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Figure  8 - X-ray  Reflection  Topographs  Correlated  to 
Rocking-Curve  Profiles  in  Figure  6.  First  image  corresponds 
to  rotation  over  entire  reflecting  range;  topographs  num- 
bered 1-10  and  1-8  correspond  to  angular  positions  indicated 
on  profiles  6(a)  and  6(b),  respectively. 
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Moreover,  the  individual  substructural  domains  in  the  surface 
layer  remained  in  reflecting  position  over  a much  larger  angular 
range  than  those  in  depth.  For  example,  the  large  reflecting 
domain  marked  by  an  arrow  went  out  of  reflection  after  an  angular 
rotation  of  only  360  seconds. 

The  structural  investigation  of  the  fatigued  polycrystalline 
2024  alloy  was  carried  out  by  using  a special  X-ray  double 
crystal  di^fr^ctygneter  technique  development  for  polycrystalline 
materials.  1 ' 3~  Item  (a)  of  figure  9 shows  the  arrays  of  spots 
corresponding  to  the  different  grain  reflections  of  the  alloy 
prior  to  fatigue  cycling. 
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Figure  9 - Rocking  Curves  of  Grain  Reflections  of  2024  A1 
prior  to  Cycling.  Cu  K monochromated  radiation. 
Angular  specimen  settings  5*  of  arc. 

(a)  Overall  survey  of  grain  reflections 
(b)  Detail  of  range  of  reflection  of  grain  A. 


The  discrete  angular  rotation  of  the  specimen  between  consecutive 
spot  reflections  was  5 minutes  of  arc.  Because  a crystal  mono- 
chromated beam  impinged  on  the  specimen  having  an  average  grain 
size  of  about  100  dm,  relatively  few  grains  were  in  reflecting 
position.  Consequently,  the  spot  reflections  along  the  Debye- 
Scherrer  ring  appeared  distinct  and  well  separated.  The  rocking 
curves  at  high  azimuthal  angles,  that  is,  at  locations  far 
removed  from  the  equatorial  reflections,  became  artificially 
broadened  due  to  tfrf  contribution  of  the  vertical  divergence  of 
the  incident  beam,  which  is  unaltered  by  reflection  from  the 
first  crystal,  it  is  most  convenient,  therefore,  to  focus 
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attention  on  the  equatorial  reflections,  such  as  spot  reflection 
A in  figure  9*  since  for  these  the  contribution  of  the  vertical 
divergence  to  broadening  is  zero.  It  can  be  seen  that  for  grain 
A the  angular  range  of  reflection  was  traversed  after  five  spot 
settings,  thus  amounting  to  about  25  minutes  of  arc  (item  (b)  of 
figure  9).  Arrays  of  spots  pertaining  to  other  grain  reflections 
exhibited  similar  or  even  smaller  values.  By  contrast,  the  grain 
reflections  in  item  (a)  of  figure  10  pertaining  to  the  cycled 
specimen  not  only  exhibited  a much  larger  angular  range,  mani- 
fested by  the  extended  arrays  of  spot  reflections,  but  also 
showed  the  typical  image  distortions  characteristic  of  a sub- 
structure induced  by  deformation.  These  two  aspects  can  be  seen 
from  the  equatorial  grain  reflection  B (item  (b)  of  figure  10 J 
which,  like  the  grain  reflection  A,  pertains  to  the  (311)  reflec- 
tion. 
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Figure  10  - Rocking  Curves  of  Grain  Reflections 
of  Cycled  2024  A1  (N  = 104 ) . 

Cu  K monochromated  radiation. 

Angular  specimen  setting  5'  of  arc. 

(a)  Overall  survey  of  grain  reflections. 

(b)  Detail  of  range  of  reflection  of  grain  B. 


Grain  B exhibited  two  multipeaked  arrays  displaced  with  respect 
to  each  other.  The  total  halfwidth  of  the  rockina  curves  cov- 
ered about  18  spot  settings  or  18  x 5 minutes  or  ~ 1.5  degrees  of 
arc.  The  multipeaked  intensity  distribution  was  due  to  the  sub- 
structural  lattice  comains  which  were  sequentially  brought  into 
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and  out  of  reflecting  position  as  the  grains  were  rotated  through 
the  angular  range  of  reflection.  The  multipeaked  nature  of  the 
intensity  curve  was  quite  analogous  to  that  obtained  for  large 
tensile -deformed  aluminum  single  crystals  shown  in  figure  6. 

After  removal  of  a surface  layer  of  100  Urn  from  the  cycled  spec- 
imen the  image  distortions  of  the  spot  reflections  decreased  and 
the  average  halfwidth  value  of  the  rocking  curves  £ was  diminished 
as  may  be  seen  from  figure  11  (viz,  grain  reflection  C). 
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Figure  ll  - Rocking  Curves  of  Grain  Reflections  of 
Cycled  2024  A1  after  Removal  of  100  Mm  Surface  Layer. 
Cu  K monochromated  radiation. 

Angular  specimen  settings  5*  of  arc. 

(a)  Overall  survey  of  grain  reflections. 

(b)  Detail  of  range  of  reflection  of  grain  C. 


Typical  values  were  between  40  and  50  minutes  of  arc.  Using  the 
average  3 values  for  the  cycled  specimen  of  90  minutes  at  the  sur- 
face and  45  minutes  at  a depth  of  100  Mm,  one  obtains  a decrease 
in  the  density  of  excess  dislocations  from  23.1  x 109  to 
5.8  x 109  (i.e.,  a reduction  by  a factor  of  about  4).  An  examina- 
tion of  the  specimen  after  removal  of  200  Mm  gave  an  excess  dis- 
location density  of  the  same  value  as  that  at  100  Mm.  The  excess 
dislocation  density  as  a function  of  depth  is  given  in  figure  12. 
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DISTANCE  FROM  SURFACE 


Figure  12  - Excess  Dislocation  Density 
of  2024  Al  Specimen  Fatigued  104  Cycles 
at  0.5  Yield  Strength 

MAGNETIC  PERMEABILITY 

Representative  oscilloscope  traces  obtained  by  the  first 
method  described  above  for  a 4130  steel  specimen  that  had  been 
strained  0.4$  is  shown  in  figure  13-  The  amplitudes  A and  the 
phases  of  the  induced  voltage  $ are  given  in  table  2.  These 
data  show  that  after  straining  0.4$  the  amplitude  of  the  induced 
voltage  A decreased.  However,  the  phase  angle  $ between  the 
input  and  output  voltage  remained  relatively  unchanged.  These 
measurements  were  found  to  be  sensitive  to  the  relative  disposi- 
tion of  and  Lg  (see  figure  2).  Although  care  was  taken  to 
assure  that  the  contribution  of  this  extraneous  factor  was  mini- 
mal, it  was  decided  to  verify  the  findings  up  to  this  point  and 
to  attain  a potentially  more  discriminating  technique  for  observ- 
ing phase-angle  changes.  Consequently,  a specimen,  strained 
plastically  approximately  2$,  was  measured  both  prior  to  and  sub- 
sequent to  deformation  by  using  the  impedance  measuring  arrange- 
ment shown  in  figure  3.  The  resultant  uncorrected  data  are  shown 
in  table  3«  Tedious  but  straightforward  corrections  for  the 
resistance  of  the  inductors  at  the  various  frequencies  as  well  as 
the  stray  inductance  of  the  standard  resistors  must  be  made.  The 
basic  data  for  these  corrections  are  summarized  in  table  4.  It 
is  evident  that,  confirming  the  earlier  findings  for  the  specimen 
strained  0.4$,  the  deformation  resulted  in  a marked  decrease 
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in  measured  inductance  at  very  low  frequencies  coupled  with  a 
decrease  in  the  measured  resistance  at  the  higher  frequencies. 


Figure  13  - (a)  (b)  - Trace  of  Voltage  across  Lp 
versus  Current  through  Li; 

(c)  ( d)  - Trace  of  Integral  of  Voltage 
across  Lp  versus  Current  through  Li 


( 
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TABLE  2 - AMPLITUDE  AND  PHASE  OF 
INDUCED  VOLTAGE  FOR  UNSTRESSED  AND  STRESSED 
SPECIMENS  OF  4130  STEEL  STRAINED  0.4$ 


A(l)* 

mv 

FjTftai 

EH 

1 

5 

6 

4.5 

0.81 

0.83 

10 

11 

9-4 

0.82 

0.84 

20 

21 

18.5 

0.75 

0.79 

50 

50 

41.5 

0.60 

0.60 

100 

76 

65 

0.405 

0.44 

200 

100 

85 

0.32 

0.32 

500 

137 

ill 

0.229 

0.214 

1000 

162 

127 

0.167 

0.157 

2000 

182 

140 

0.16 

0.128 

5000 

200 

158 

0.112 

10,000 

231 

182 

0.143 

* Peak-to-peak  amplitude  of  induced  voltage  for  a 
primary  current  of  25  mA  peak  to  peak; 

A( 1 ) unstressed,  A(2)  stressed. 

**Ratio  of  minor  to  manor  axis  of  ellipse;  convert 
to  phase  angle;  from  theory  in  appendix  C; 

0(1)  unstressed,  0(2)  stressed. 


TABLE  3 - IMPEDANCE  MEASUREMENTS  OF  4130  STEEL 
UNSTRAINED  AND  STRAINED  2$ 


8 

L,  henries 

R, 

ohms* 

Unstrained 

Strained  2$ 

Unstrained 

Strained  2$ 

5 

1.065 

0.869 

86 

100 

10 

1.015 

0.868 

92 

102 

20 

1.017 

0.861 

102 

110 

50 

0.922 

0.809 

165 

154 

100 

0.771 

0.698 

285 

258 

200 

0.624 

0.558 

455 

419 

500 

0.468 

0.412 

640 

753 

1,000 

0.377 

0.329 

1,350 

1.190 

2,000 

0.306 

0.266 

2,080 

1.828 

5,000 

O.227 

0.200 

4,000 

3,400 

10,000 

0.186 

■ ...  — — — ■ 

0.166 

6,900 

5,800 

♦Readings  must  be  corrected  for  resist  nee  of  standard  inductors, 
also  one  must  subtract  fixed  resistor  of  105  ohms  in  series 
with  Li  inserted  to  allow  balance  at  lower  frequencies. 
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TABLE  4 - CORRECTION  FACTORS  FOR  SPECIMEN 


L Setting, 
henries 
x 10-5 

R 

ohms 

L Setting, 
henries 
x 10-2 

R 

ohms 

L Setting, 
henries 
x 10-1 

R 

ohms 

1 

0.26 

1 

0.92 

1 

8.55 

2 

0.55 

2 

1.54 

2 

15.1 

5 

0.54 

5 

2.58 

5 

25.6 

4 

0.48 

4 

2.80 

4 

28.4 

5 

0.68 

5 

5.67 

5 

56.9 

6 

0.78 

6 

4.50 

6 

45.5 

7 

0.96 

7 

5-15 

7 

52.0 

0 

1.02 

8 

5-97 

8 

61.6 

9 

1.09 

9 

6.57 

9 

68.2 

10 

1.26 

10 

7.42 

10 

76.6 

Inductance  of  Standard  Resistors  (at  10  kHz’) 

R Setting 
ohms 

Inductance 

mH 

R Setting 
ohms 

Inductance 

mH 

100 

0.5 

900 

2-9 

200 

0.7 

1,000 

5.5 

500 

1.0 

2,000 

6.5 

400 

1.5 

3,000 

9-5 

500 

1.6 

4,000 

13.0 

600 

1.9 

5,000 

16.0 

700 

2.2 

6,000 

19.0 

800 

2.5 

7,000 

22.0 

To  determine  the  effect  on  the  impedance  of  the  reduction  in 
diameter  of  the  specimen  strained  2$,  an  unstrained  specimen  was 
electropolished  co  reduce  the  diameter  the  same  amount  (0.25  mm). 
Impedance  measurements  were  made  before  and  after  the  metal 
removal.  From  these  data  it  was  apparent  that  the  r cc_ jticr.  in 
the  diameter  did  not  account  for  the  major  changes  observed  in 
the  strained  ( 2 %)  specimens. 

ACOUSTIC  EMISSION 

The  acoustic  emission  behavior  was  investigated  on  speci- 
mens of  4150  steel  that  had  been  fatigued  at  90  and  125  ksi 
after  1 cycle  and  after  25$,  50$,  and  75$  of  the  fatigue  life. 

It  was  quite  apparent  in  all  cases  that  emissions  occurred  at 
stresses  quite  far  below  those  at  which  the  specimens  had  been 
fatigued.  This  portion  of  the  investigation  was  not  conducted 
in  sufficient  depth  to  determine  the  reasons  for  the  nonexis- 
tance  of  the  Kaiser  effect;  however,  1.0m  the  acoustic  emission 
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research  currently  under  way  by  Carpenter,  “ it  appears  that  the 
magnetostriction  properties  of  steels  are  important.  Such  prop- 
erties cause  the  domain  configuration  to  be  determined  in  part  by 
the  internal  stress  fields.  Applying  a load  will  cause  rapid 
changes  in  the  domain  configuration  in  the  manner  of  the  Bark- 
hausen  effect.13  These  changes  will  be  the  source  of  elastic 
waves  via  the  changes  in  the  magnetic  structure  that  will  be 
sensed  by  the  pickup  transducer.  Still  another  possibility  for 
the  low  stress  emissions  is  relaxation.  However,  it  appears  that 
the  emissions  occur  at  too  low  a stress  relative  to  the  prestress 
to  be  accounted  for  by  such  as  mechanism. 

DISCUSSION 

The  X-ray  line  broadening  and  topography  study  shows  that 
tensile  deformation  in  such  diverse  materials  as  silicon  and 
aluminum  resulted  in  a preferential  deformation  of  surface  layers 
compared  to  that  of  the  bulk  material.  in  aluminum,  the  deformed 
surface  layer  exhibited  not  only  smaller  misaligned  lattice 
domains  than  those  in  depth  but  also  the  larger  angular  range 
over  which  the  individual  substructural  reflections  persisted  was 
indicative  of  a greater  density  of  excess  dislocations  within 
these  domains.  The  substructural  features  of  the  surface  layer 
may  be  likened  to  a lattice  "debris,"  whereas  those  in  depth  are 
more  like  the  substructure  of  the  annealed  material.  The  differ- 
entiation is  attributed  to  preferred  activation  of  dislocation 
sources  located  near  the  surface8 a and  to  intense  interaction  of 
moving  dislocations  on  multiple  slip  systems.  Kitaj ima3 0 31  c?~3 
to  similar  conclusions  using  etch-pit  techniques  for  the  charac- 
terization of  the  dislocation  structure  in  deformed  copper. 

Kramer  9 '3  ~a 4 has  reported  previously  that  work  hardening  does 
not  occur  uniformly  throughout  a uniaxially  strained  specimen  and 
is  greater  in  the  surface  layer  than  in  the  interior.  He  reporte 
that  this  surface  layer  also  forms  in  stage  I of  fee  metals  where 
only  single  slip  is  supposed  to  operate. 

The  preliminary  results  obtained  from  the  fatigue  cycled 
A1  2024  alloy  indicate  a strong  adherence  to  the  "debris"  con- 
cept.39 The  X-ray  images  derived  from  the  reflecting  grains  of 
the  undeformed  material  were  uniform  and  passed  quickly  in  and 
out  of  reflection,  viz,  20  to  25  minutes  of  arc.  In  contrast, 
when  a specimen  was  subjected  to  10,000  fatigue  cycles  at  a 
stress  amplitude  amounting  to  0*5  of  the  yield  stress  of  the 
material,  the  images  persisted  over  a large  angular  range,  viz, 

90  minutes  of  arc,  and  were  broken  up  indicating  the  formation  of 
a highly  misoriented  surface  layer.  In  addition,  the  intensity 
of  the  spots  increased  significantly  over  their  undeformed  coun- 
ter parts  signifying  an  increase  in  integrated  intensity  with 
the  broadening  of  the  associated  rocking  curves.  The  reflec- 
tions emanating  from  greater  depths  of  the  material  and  identi- 
fied by  their  weakness  relative  to  surface  reflections  of  the 
same  crystallographic  order  were  found  to  be  highly  regular  and 
short  in  reflecting  range.  This  suggests  that  the  bulk  material 
did  not  incur  as  much  damage  as  the  surface  layer  as  a result 
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of  fatigue  cycling.  The  X-ray  results  obtained  after  a surface 
layer  of  100  Mm  had  been  removed  confirmed  the  greater  degree  of 
crystal  perfection  since  the  reflection  curves  became  narrower, 
vis,  45  minutes  of  arc,  and  the  break-up  of  the  grain  reflections 
into  substructural  units  declined.  Electron  microscope  studies 
by  Duquette  and  Swann  provided  similar  evidence  of  high  surface 
layer  dislocation  densities  in  a low-strain,  fatigue-cycled 
aluminum  alloy.'1' 

With  reference  to  the  cause  of  the  magnetic  measurement 
changes  upon  deformation  of  specimens  strained  0.4$  and  2$,  the 
data  in  table  3 show  that  there  is  a pronounced  decrease  in  the 
measured  inductance  at  the  low  frequency  of  5 Hz  after  deforma- 
tion. From  equation  (l4a)  and  figure  1 it  may  be  seen  that 
t^eff.r  rs  1 for  very  low  frequencies  (small  f/fg);  consequently, 
a decrease  in  inductance  must  be  attributed  to  a decrease  in  M. 

A more  quantitative  estimate  of  f/fg  for  the  specimens  in  ques- 
tion is  contained  in  appendix  E.  Tnese  estimates  confirm  that 
f/f„  is  very  small  for  f = 5 Hz.  Some  supplemental  measurements 
of  Specimen  resistance  were  made  by  using  a Kelvin  bridge.  The 
resistance  was  found  to  be  approximately  4.3  milliohms  for  all 
specimens  and  was  relatively  unchanged  by  plastic  deformation, 
thereby  confirming  the  conclusion  that  measurement  changes  were 
due  to  changes  in  the  initial  permeability  M0.  The  changes  in 
permeability  close  to  the  surface  may  be  estimated  with  the  aid 
of  equation  (14b)  applied  to  the  highest  frequency  measurements, 
equation  (14a)  being  insensitive  to  higher  frequencies  results 
due  to  the  dominance  of  the  1 - h term.  From  this  equation  it  is 
easily  seen  that 

AR  ’ /R ' = AM/M  t AMeff.^Meff.i  (15) 

where  R’  = R - R . From  equation  (13)  and  table  1, 

AMeff # i/^eff . i may  be  related  to  Am/M.  At  the  very  highest  fre- 
quencies AMeff.i  is  independent  of  M.  At  lower  frequencies  the 
relationship  may  be  approximated  as  is  shown  in  appendix  F.  The 
main  point  is  that  a decrease  in  M at  these  frequencies  is  par- 
tially compensated  by  an  increase  in  Meff.i  so  that  the  overall 
change  in  R'  is  diminished.  By  applying  this  to  the  data  the 
change  in  M may  be  estimated  and  compared  to  the  change  in  M at 
5 Hz  (equation  (14a)  as  developed  in  appendix  G).  The  results 
of  these  calculations  are  summarized  in  table  5*  From  these  it 
is  seen  that  the  plastically  deformed  specimens  seen  to  exhibit 
substantially  larger  changes  in  M at  higher  frequencies. 
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TABLE  5 

CHANGE  IN  U at 
LOW  AND  HIGH  FREQUENCIES 


Strain  % 

Frequency 

Ab/U 

0.4 

5 Hz 

-0.27 

10  kHz 

-0.42 

2 

5 Hz 

-0.20 

10  kHz 

-0.35 

Some  mention  should  be  made  of  the  effective  depth  of  perme- 
ability measurement.  Calculations  of  skin  depth  as  a function  of 
f/fg  for  specimens  of  cylindrical  geometry  have  appeared  in  the 
literature.16  Applying  the  results  to  the  specimens  reported 
herein  yields  the  skin  depths  listed  in  table  1-F  of  appendix  F. 
It  is  seen  that  the  skin  depth  for  the  highest  frequency  70kHz 
is  approximately  125  Wn#  which  depth  would  appear  to  be  the 
desired  order.  However,  it  should  be  realized  that  depth  varies 
inversely  as  pl/2.  Consequently,  a decreasing  p will  result  in 
a greater  skin  depth  and  the  consequent  change  in  measured  param- 
eters counter  to  that  produced  by  the  decrease  in  U.  This  makes 
the  overall  technique  somewhat  less  sensitive  than  it  would 
appear  to  be  at  first  sight.  However,  the  limited  data  obtained 
thus  far  would  indicate  that  radial  variations  of  U are  obser 
vable . 

Not  much  mention  has  been  made  of  coercivity  measurements, 
although  the  oscilloscope  displays  shown  in  figure  13  would 
allow  its  estimation  throughout  the  entire  cross  section.  At 
higher  frequencies,  eddy  currents  confounded  the  display  to  a 
point  where  simple  interpretation  was  not  possible.  Further- 
more, the  effects  of  the  demagnetizing  field  at  the  ends  of  the 
specimen  are  not  easily  ascertained,  although  it  may  be  shown 
shown1  ,'i7  that  the  effects  on  measured  permeability  were  negli- 
gible for  the  specimen  geometries  and  ranges  of  permeability 
under  consideration.  Consequently,  this  particular  measurement 
technique  offered  little  in  the  way  of  measurement  of  th^  coer- 
civity in  the  surface  layer  and  so  was  not  pursued  further. 

CONCLUSIONS 

• From  the  X-ray  line  broadening  and  B-B  topographic 
measurement  it  was  quite  evident  that  the  excess  dislocation 
density  in  the  surface  layer  was  much  greater  than  that  in  the 
bulk  for  specimens  that  had  been  fatigued  or  strained  unidirec- 
tionally . 


• From  the  measurements  on  plastically  deformed 
steels  that  showed  that  the  initial  permeability  decreased  at 
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high  frequencies  (10  kHz),  it  could  also  be  inferred  that  the 
dislocation  density  was  larger  in  the  surface  layer  than  in  the 
bulk. 


• Because  emissions  occurred  at  relatively  low 
stresses  after  fatiguing  at  high  stress  amplitudes,  it  appears 
that  acoustic  emission  cannot  be  used  to  detect  fatigue  change  in 
ferromagnetic  materials.  While  such  a technique  may  possibly  be 
used  in  the  laboratory,  it  could  not  be  used  in  the  field  unless 
the  magnetic  domains  can  be  prevented  from  growing.  There  is  a 
possibility  that  acoustic  emission  may  be  useful  for  nonmagnetic 
materials  - provided  other  emission  sources  do  not  operate. 

FUTURE  WORK 

During  the  following  year  emphasis  will  be  placed  on  the 
X-ray  rocking  curve  of  polycrystalline  aluminum  and  magnetic 
permeability  of  steel.  The  line  widths  will  be'  measured  as  a 
function  of  the  fraction  of  fatigue  life.  For  the  magnetic 
permeability  portion  of  this  effort  measurement  at  higher 
frequency  (>10  kHz)  will  be  investigated  to  obtain  improved 
sensitivity. 
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APPENDIX  A 


RELATIONSHIP  OF  INTERNAL  STRESS  TO  INITIAL  PERMEABILITY 


I 
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Van  Nostrand 


This  parallels  the  simplified  derivation  given  in  the  early 
literature  (reference  (a)).  First  consider  a magentic  domain 
over  which  there  is  a uniform  internal  stress  in  the 

absence  of  an  applied  field,  the  domain  magnetization  will 
assume  a minimum  energy  direction,  which  direction  is  parallel 
to  t ^ for  materials  with  positive  magnetostriction.  If  a mag- 
netic field  H is  applied  an  an  angle  90  with  respect  to  t^, 
then  the  magnetization  of  the  domain  Is  will  rotate  partially 
into  this  direction  to  some  equilibrium  angle  0 ( see  figure  1-A) . 


Among  H,  Is,  and  T 


The  angle  6 is  that  which  results  in  the  total  domain  energy 
being  a minimum.  It  can  be  shown  (reference  (b))  that  the  total 
ebergt  E has  two  components,  the  strain  energy  E<ri  and  the 
energy  of  a dipole  in  a magnetic  field  Eh#  and  that  these  are 
given  by 


I 


ETi  = 3/2  (\g  TiSin20) 


(A-l) 


I 
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and 


eh  = -Hiscos  (e0  - e)  (A-2) 

where  Xa  is  the  magnetostrictive  coefficient  at  saturation  and 
Is  is  the  saturation  magnetization  per  unit  volume.  The  condi- 
tion for  equilibrium  is  that 


dE/d6  = 0;  where  E = ETi  + EH  (A-3) 

Applying  this  condition  results  in 

3XS  i sin  9 cos  6 = HIS  sin  (0O  - 9)  (A-4) 

For  small  H and  consequent  small  9,  sin9  = 9,  cos9  ~ 1 and 
H0  ~ 0,  so  that  equation  (A-4)  becomes 

0 = [HIS  sin0oj/6xsxi  (A-5) 


To  calculate  the  permeability,  we  first  note  that  for  small 
applied  fields  the  magnetization  I in  the  direction  of  the 
applied  field  H is  given  by 


I = kH 


(A— 6 ) 


where  k is  the  magnetic  susceptibility.  Since 

B = H + I (A-7) 

by  definition  and  since  the  permeability  M.  = B/H  by  definition 
as  well,  it  is  seen  from  equations  (A-6)  and  (A-7)  that 

M - 1 = k = dl/dH  = (dl/d0 ) x ( d 0/dH ) (A-8) 
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Since  the  component  of  l„  in  the  direction  of  the  field  is  given 

S 

by 


for  an  internal  stress  at  an  angle  6 with  respect  to  the  applied 
field  H.  If  we  average  this  expression  over  all  possible  values 
of  0O  from  0°  to  180°  we  obtain  sin2  0O  = 1/2  over  this  interval, 
so  that 

yx  - 1 = ig^XgTi  (A-13, 


While  this  is  only  a very  approximate  expression,  it  tells  us 
that  permeability  varies  inversely  with  internal  stress  and  that 
absence  of  stress  should  result  in  very  large  permeabilities. 
This  does  occur  for  highly  oriented  microstructures  wherein  me 
grains  are  oriented  along  the  "easy"  axis  of  magnetization.  In 
a randomly  oriented  metal  the  magnetizing  field  must  overcome 
the  restoring  forces  of  crystalline  anisotropy  as  well  so  that, 
even  in  the  absence  of  internal  stress,  the  permeability  would 
not  be  as  high  as  predicted  by  equation  (A-13)»  However,  it 
shows  that  the  presence  of  internal  stress  should  decrease  the 
permeability. 


APPENDIX  B 

PENETRATION  OF  A VARYING  MAGNETIC  FIELD  INTO  A CONDUCTOR 
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This  derivation  is  given  primarily  to  expose  the  physical 
assumptions  underlying  most  of  the  statements  regarding  the  so 
called  "skin  effect"  and  thereby  give  the  understanding  required 
to  appreciate  the  results  of  the  investigations  reported  herein. 
We  start  with  Maxwell's  equations  for  a homogeneous,  isotropic, 
linear  medium 


V x E 


V x H 


dB 
~ dt 

(B-l) 

dD 

= dt  + j 

(B-2) 

where  E is  the  electric  field  intensity,  H the  magnetic  field 
intensity,  D the  electric  displacement,  B the  magnetic  induction, 
j the  current  density,  and  t the  time.  Units  need  not  be  of  con- 
cern at  the  moment.  The  major  assumption  made  at  this  point  is 
that  for  the  electrical  conductivities  and  frequencies  of  inter- ;t 
we  may  ignore  the  displacement  current  in  comparison  to  the 
conduction  current  j ? or  dt 

V x H = j (B-3) 

In  addition,  it  is  assumed  that  Ohm's  law  holds  for  the  medium? 
i.e. , 

j = OE  (B-4) 

where  o is  the  electrical  conductivity,  and  further  that 

B = PH  (B-5) 


where  p is  the  magnetic  permeability.  (This  holds  even  for 
magnetized  ferromagnetic  materials  provided  H is  sufficiently  small, 
and  the  correct  value  is  used  for  p.)  With  these  assumptions, 
equations  (B-l)  and  (B-2)  become 


7 X E 


PdH 

dt 


(B-6) 


7 x H = OE 


(B-7) 
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B-l 


Taking  the  curl  of  both  sides  of  both  equations  and  using  the 
vector  analysis  identity 


V X V x = VV.  _V»  (B-8) 

and  by  manipulating  to  obtain  isolated  equations  for  E and  H the 
following  partial  differential  equation  is  obtained  for  H 

7aH  = n<J  (B-9) 


This  is  the  vector  equivalent  of  the  heat  conduction  or  diffusion 
equation  which  already  has  been  solved  for  many  practical  cases 
of  interest.  With  the  aid  of  figure  1 of  the  text  the  simplest, 
planar  solution  may  be  derived.  H is  assumed  to  be  parallel  to 
the  surface  of  a semi-infinite  metallic  region  and  to  have  the 
following  simple  time  dependence  external  to  the  metal  with  no 
external  spatial  dependence; 

Hz(x)  = Hqcos  ( 2n  ft)  ; x>o  (B-10) 

Direct  substitution  in  equation  (B-9)  for  which  we  also  assume  Hx 
f Hy  = 0 verifies  that  the  following  expression 

Hz  = H sx/s  cos  (2  nft  - - \ (B-ll) 

S' 

is  a solution  in  the  planar  case  satisfying  the  boundary  condi- 
tion (reference  (17),  provided  that 

s = (TTpof)*15  ( B— 12) 

s is  obviously  the  depth  at  which  H decreases  to  1/e  times  its 
value  at  the  surface,  x = 0. 


In  the  cylindrical  case,  the  boundary  equation  (B-10)  remains 
the  same  for  r>r Q,  the  radius  of  the  cylinder.  Equation  (B-9)  in 
cylindrical  coordinates  for  Hr  = H = 0 becomes 


13  A to 


r dr\ 


?> 


HO 


at 


(B-13) 


From  the  standard  approach  of  separation  of  variables  we  assume 


Hjj  = R(r)T(t) 


which  yields 


3_  dR\  no  dT 
rR  dr  \ dry  T dt 


(constant) 


(B-14) 
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The  equation  in  R is,  upon  appropriate  transformation  of  variables, 
the  zeroth  order  Bessel  equation.  Consequently,  the  solution  for 
the  spatial  dependence  of  H turns  out  to  be 


H 


Ho 


bera 

bera 


(2  6r/d)+  bei 
9 +beia  (0) 


( B— 1 5) 


where  the  ber  and  bei  functions  are  complex  combinations  of  ' 
Bessel's  functions  and 


./  rrpifO 

0 = d V -=2~  (B-16) 

where  d is  the  diameter  of  the  cylindrical  specimen. 

Equation  (B-15)  provides  the  basis  for  computing  the  effec- 
tive permeability  of  a cylindrical  specimen.  The  effective 

permeability  is  defined  as  that  value  that  yields  the  total  flux 
cp  through  the  cross  section  of  a cylindrical  specimen  in  the 
following  equation 

cp  = nrca  P Meff  Ho?  rc  = d/2  (B-17) 

However,  following  the  standard  derivation  (reference  (a)),  the 
total  flux  is  obtained  by  summing  the  flux  through  each  cylindrical 
shell  of  thickness  dr  as  expressed  by  the  equation  below 


r H dr 


(B-18) 


where  H is  given  by  equation  (B-15) . From  equations  (B-17)  and 
(B-18)  we  have 


eff 


H. 


L-  fr‘ 

>ro  J 


rHdr 


(B-19) 


Upon  integration,  equation  (B-19)  yields  (reference  (b) ) 

_ 2 


eff 


J1(a9/2) 


Krc  Jo(a0/2) 


(B-20) 


where  JQ  and  are  Bessel  functions  of  the  zeroth  and  first 
order  and  a as  a complex  quantity  is  given  by 
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a = fT.  where  j = V~-l 


(B-21) 


It  follows  that  Meff  is  a complex  quantity.  In  physical  terms, 
this  means  that  the  flux  through  the  specimen  is  out  of  phase 
with  the  applied  magnetizing  field. 


APPENDIX  C 

COMPUTATION  OF  PHASE  FROM  AN  ELLIPTICAL 
LISSAJOU  FIGURE 

Assume  that  the  x and  y axis  deflections  of  an  oscilloscope 
display  are  described  by  the  following  time  parametric  equations 

x = rsin  (tot  + a-^)  (C-l) 

y = rsin  (urt  + a2)  (C-2) 

That  is,  each  axis  deflection  is  sinusoidal,  of  equal  amplitude, 
and  of  differing  phase.  Upon  some  manipulation  and  the  use  of  a 
standard  trigonometric  identity,  these  become 

x/r  = sinwtcosaq  + coswtsina^  (C-3) 

y/r  = sinwtcosa.2  + coswtsina2  (C-4) 

The  following  manipulations  are  then  required  to  eliminate  time 
from  these  two  equations : 

1.  Multiply  equation  (C-3)  by  sina.2  and  equation  (C-4) 

by  sinaq. 

2.  Subtract  the  modified  equation  (C-4)  from  the  modi- 
fied equation  (C-3) . 

3.  Square  the  resultant  equation. 

4.  Multiply  equation  (C-3)  by  cosa2  and  equation  (C-4, 

by  cosaq. 

5.  Subtract  the  modified  equation  (C-4)  from  modified 
equation  (C-3) . 

6.  Square  the  resultant  equation. 

7.  Add  the  equation  obtained  by  step  3 to  the  equation 
obtained  by  step  6. 

By  judicious  use  of  the  trigonometric  identity 

sinae  + cosa6  = , 1 (C-5) 

the  equation  resulting  from  step  7 becomes 


I 
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j 2 

x + y 2xycos  (a2  -a^)  ^ 

| r8  sin8 ( a2  -ai) 

(C-6) 

This  equation  may  be  shown  to  be  that  of  an  ellipse  with  the  major 
axis  inclined  45°  to  the  x axis.  Starting  with  the  equation  for 
an  ellipse  with  the  major  axis  along  the  x axis 

x8/aa  + y8/b8  = 1 

(C-7) 

where  a and  b are  the  lengths  of  the  major  and  minor  axes, 
respectively.  The  following  rotational  transformation  will  generate 
the  equation  of  a 45°  inclined  ellipse 

x ' = x - y 
y'  = x + y 

(C-8) 

resulting  in  (dropping  the  prime  sign) 

(x+y)  a/as  + (x-y) 3/b3  = 2 

(C-9) 

as  the  equation  for  a 45°  inclined  ellipse.  Comparing  this  to 
equation  (C-6)  it  may  easily  be  shown  that 

a8  = (sin8  (a2  -a1)r8)/(l  - cos(or2  - <*i)) 

(C-10) 

b8  = (sin8  (a2  -ai)r8)/(l  + cos(a*2  ~ al)) 

(C-ll) 

so  that,  dividing  equation  (C-ll)  by  equation  (C-10)  and 
for  cos  (a2  -ai) , we  obtain 

solving 

cos  (a2  -al)  = ( 1-  (b/ a)  8 ) /(  1+  (b/ a)  8 ) 

(C-12) 

Thus,  to  measure  the  phase  difference  between  the  x and 
tions  where  they  are  of  the  same  amplitude,  measure  the 
majo’"  axes  of  the  generated  ellipse  and  substitute  into 
(C-12) . 

y deflec- 
minor  and 
equation 
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APPENDIX  D 

TRACE  OF  B/H  CHARACTERISTIC 

The  emf  generated  in  of  figure  2 of  the  text  is,  by 
Faraday's  law  of  induction,  given  by 

e2(t)  = Kdcp/dt  (D-l) 

where  K is  a constant  and  cp  is  the  total  flux  of  B through  the 
coil,  i.e.. 


Hence,  if  we  subject  e2  to  integration  via  the  simple  RC  integrat- 
ing circuit  shown,  we  obtain  the  following  for  the  vertical  deflec 
tion  of  the  oscilloscope 

y ( t)  = K'cp(t)  (D-3) 

If  we  assume  that  the  primary  L-^  generates  a constant  magnetic 
field  across  L2  that  is  proportional  to  the  primary  current,  then 

x ( t)  = K"H(t)  (D-4) 

where  K"  is  another  proportionality  constant.  Eliminating  the 
time  t between  equations  (D-3)  and  (D-4)  will  result  in  the 
following  relationship  between  the  x and  y deflections 

B.dA  = G (H) 

If  B is  constant  across  the  area,  this  reduces  to 

G (H) 

B = A 

which  is  the  B/H  characteristic  of  interest  aside  from  a scale 
factor.  In  fact,  H is  not  constant  across  L2  and  for  certain 
instances,  neither  is  B,  if  material  properties  vary  with  radial 
distance  from  the  surface.  To  this  extent,  the  circuit  shown  will 
only  give  a sort  of  average  B/H  characteristic. 
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APPENDIX  E 
CALCULATION  OF  f 

9 
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By  definition  (references  (a)  and  (b) ) , f is  given  by  the 
following  expression  for  a cylindrical  specimen 

fg  = 5066/ (pod3)  Hz  (E-l) 

where  the  units  of  the  variables  are 

P = relative  permeability,  dimensionless 

o = electrical  conductivity,  meter/ohm-mm3 

d = diameter  of  cylinder,  cm 

Based  on  the  measurements  made,  the  following  values  were  used  to 
compute  fg: 

P = 47 

o = 3.5  meter s/ohm-mm3 
d = 1.2  cm 

These  result  in  a value  of  21  Hz.  The  basis  for  each  of  the  above 
will  now  be  presented. 

The  value  of  p is  approximately  the  same  as  reported  else- 
where for  similar  steels  and  is  based  on  the  change  in  inductance 
of  coil  L-^  at  low  frequencies.  The  value  of  o is  based  on  a 
Kelvin  bridge  measurement  of  specimen  resistance.  At  low 
frequencies  Peff>r  approaches  unity  so  that 

L = L0  (1-n+rip)  ( E-2) 

the  inductance  of  L^  devoid  of  a specimen  was  measured  as  being 
0.093  henry,  while  for  a multilayer  coi]  of  the  dimensions  of  L^ 
the  average  value  of  T]  was  calculated  to  be  "H  = 0.23.  Based  on 
these  and  the  measured  value  of  inductance  at  5 Hz  after  insertion 
of  the  specimen,  the  value  of  P was  calculated  from  equation  (E-2) . 


APPENDIX  F 

CALCULATION  OF  AM/M  FROM  HIGH-FREQUENCY  MEASUREMENTS 


To  relate  A(aef f . i/^ef f . i to  AM/M,  as  is  required  if  we  wish 
to  use  equation  (10)  of  the  text,  we  first  go  to  equation  (13) 
from  which  it  is  easily  seen  that 


f' 


f/fg 


MO  da  f 
5066 


(F-l) 


Consequently, 

Af'/f'  = AM/M  (F-2) 

From  a graphical  plot  of  table  1 of  the  text,  we  may  numerically 
derive  the  quantity  n defined  as 

n = (AMeff  i/Af')  x (fVMeff>i)  (F-3) 

for  each  f'  by  estimating  the  slope  at  each  f ' . From  eauation 
(F-2)  it  is  also  seen  that,  by  substituting  in  equation  (F-3) , 

n = (AMeff>i/Meff#i)/AM/M  (F-4) 

so  that 

Ahieff.i/^eff.i  = n x (AM/M)  (F-5) 

and  from  equation  (15)  of  the  text, 

r'/R'  = ( 1+n)  x (AM/M)  (F-6) 

Note  that  n is  a negative  number. 

To  complete  the  basis  for  this  analysis,  n has  been  computed 
for  the  frequencies  of  interest  from  equation  (F-3)  ana  table  1 
of  the  text.  These  are  summarized  in  table  1-F  for  the  higher 
frequencies  of  interest. 


TABLE  1-F 

TABULATION  OF  n FOR  PLASTICALLY  DEFORMED 
4130  STEEL  SPECIMENS 


Frequency 

khz 

f/£g 

^eff.i 

f,/M  __  . 
ef  f . i 

AMeff.i 

- 

Ini 

S 

mm 

Af  • 

10. 

500 

0.0664 

7,500 

7 . 2x10~b 

0.54 

0.416 

1,000 

0.0447 

21,300 

2.4x10~b 

0.51 

0.23 

2,500 

0.0421 

59,200 

4x1°-b 

0.24 

0.146 

i ns 

3.  500 

0.0387 

90. 500 

<4x10 

0.2 

0.136 

V AS 

•w*  • 
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APPENDIX  G 

LOW  FREQUENCY  CALCULATION  OF 

For  the  specimens  under  consideration,  f/fg  = 0.25  at  5 Hz. 
This  is  based  on  the  computed  value  of  fg  given  in  appendix  E. 

For  this  value  of  f/fg  it  is  seen  from  figure  1 or  table  1 of  the 
text,  that  ^eff.i  ~ 1*  80  that  equation  (14a)  of  the  text  reduces 
to 


L = Lq  (1  - ti  + tip)  (G-l) 

From  this  equation  it  is  seen  that 

AL/L  - TiAH/d-h+Tlp)  = (A*V^)/(l+(l--n)/r|p)  (G-2) 

Substituting  the  numerical  values  used  previously  for  and  n we 
obtain 


Au/u  = 1.07  (AL/L) 


(G-3) 
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